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INTRODUCTION 


Our purpose is to determine an optimal pulse shape and duration of cathodic current for 
pacing —just at its threshold- a given target element in an excitable tissue. 

The target element may be a single cell (such as nerve fiber or skeletal muscle fiber) or a 
whole network of electrically coupled cells (such as the myocardium or visceral smooth 
muscle). 

The response of the target element to the applied pulse is the all-or-none emergence of an 
action potential. 

We use the criterion of least electric energy dissipated in the tissues to select an optimal pulse 
shape amongst the different (just-threshold) cathodic ones of a given duration. 

We calculate the corresponding minimum of dissipated electric energy and then we study its 
variation with the duration of the pulse. 


METHODS 


From a phenomenological standpoint the excitation process of a given target element in an 
electrode-tissues system may be described by an “excitation functional” “©, 


For cathodic make excitation with a just-threshold pulse i(t) of duration tp, the criterion of the 
excitation functional (in the linear approximation) reduces to the following equality : 


lp 


| ¢(t.-1)-i(0)dt=4, [1] 


0 


Here qu is a certain threshold charge and g(t) is an impulsive response function studied 
elsewhere °°)-@ (Figure 1). 
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Fig. 1: The impulsive response function of an “excitation functional”. 


This function is positive between t=0 and t=tu, where tu is the useful time up to the rheobase 
(defined using strength-duration curves) ©) (Figure 2). 
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Fig. 2: The useful time up to the rheobase tu 


Furthermore, to apply [1] as a starting point, we suppose that at the beginning of the current 
pulse the excitable membranes of the target element are at their rest state. 
Applying the Cauchy-Bunyakowsky-Schwarz inequality 


feltp—t)-i(A)at <[g(r)de-[? (sae [2] 


The inequality is an equality if and only if: 


World Congress on Medical Physics and Biomedical Engineering, Sydney, 
Australia, 2003. 


i(t) = K g(tp-t) [3] 


Now from [1], for a just-threshold pulse the left-hand member of [2] must be equal to the 
constant: q;. 


Then, for a just-threshold current pulse: 


fP@a> —&— [4] 
? Je (da 


But if the biological tissues can be represented by an ohmic load of resistance Rr, the 
dissipated electric energy per pulse, can be calculated as: 


R, fi? (t) dt [5] 


RESULTS 


For a just-threshold pulse this energy will be a minimum if inequality [2] is an equality, that 
is if i(t) verifies equation [3]. Substituting [3] into [1], we determine the constant of 
proportionality K. After that, from [3] we obtain the optimal pulse shape, for cathodic make 
excitation, and for a given duration tp between 0 and tu: 


i) = _—*—.- g(,-1) [6] 


Je (t) dt 


From [4] and [5] it follows that the energy dissipated by an optimal pulse of duration tp is 
given by 


R, qa 
fe (t) dt 
0 


To use formula [6] for the optimal pulse shape and the expression [7] for the minimum of 
electric energy dissipated in 


[7] 


the tissues, both qu and g(t) must be known. In principle, they can be determined from the 
measurement of the corresponding strength-duration curves for each target element in its 
electrode-tissue system. To see how this can be done, let us consider a rectangular pulse of 
cathodic current, of amplitude io and duration tp. From the criterion given by equation [1], the 
pulse will be just threshold if and only if: 
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<4) [8] 


In this formula, tp is between 0 and ty. Now, if io=io(tp) is known (measuring the 
corresponding strength-duration curve), from equation [8] it follows that 


qd, = Jim (¢, -ip(tp)) [9] 
_ 4} 4% 
Nae Fal md 


for each tp between 0 and ty. 


In practice equation [10] shouldn’t be used to calculate g(t) form experimental data. 

The corresponding numerical problem has a very bad condition number, so that there is too 
much uncertainty in the estimation of g(t). 

To bypass this problem, we can use a closed form analytical formula with a few adjustable 
parameters to represent g(t), substitute it in formula [8] and adjust the parameters to an 


experimental strength-duration curve. 

The required formulae can be obtained from suitable mathematical models of threshold 
dynamics for fibers © or networks of electrically coupled cells © °) stimulated by external 
electrodes. 

Also, it is possible to directly introduce the linear version of the excitation functional for a 
cathodic pulse (Equation [1]) “©, and then to use a descriptive mathematical model for g(t), 
with suitable free parameters. 


DISCUSSION 


From [6] the optimal pulse shape of duration tp can be determined. 
2 
As g(t) is decreasing (between 0 and tm , see Figure 1) and convex (because ae g(t) >0), 


the optimal pulse shape is increasing and convex. 

From [7] it follows that as tp increases, the minimum dissipated energy per pulse decreases at 
its turn and reaches its minimum when tp=tu . 

The experiments show that the parameters of the strength-duration curves, and therefore qu 
and g(t), depend not only on the excitability properties of the membranes of the target 
element, but also depend on certain geometric and bioelectric parameters of the system 
formed by the electrodes, the tissues, and the target element within the tissues. 
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For example, if the target element is a fiber, they depend on the radius of the fiber. Besides, 
they depend on the distance between the cathode and the target element, and of the size and 
shape of the active electrode if it is located close enough to the target element. 

Furthermore, the dissipated energy per pulse is proportional to Rr. 

But Rr usually depends on the size and shape of the cathode (we suppose monopolar 
stimulation) and it is very sensitive to conductivity variations in the tissues located near the 
active electrode (such as the ones related with inflammatory reactions to implant electrodes). 
Therefore, the optimal pulse shape for a given duration can be modified, in the same 
electrode-tissues-target element system, due to the variation of the parameters of the system. 
A new digital stimulator for threshold research and synthesis of optimal pulse shapes was 
developed, considering these results ©), 


CONCLUSIONS 


(a) For a given duration, and a just-threshold pulse, the optimal shape is given by equation 


[6]. 


(b) The energy dissipated by an optimal pulse shape decreases as the pulse duration 
increases. 


(c) Besides being increasing and convex, there doesn’t seem to be other general properties of 
the optimal pulse shape. For each target element in its corresponding electrode-tissues 
system, the optimal pulse shape should be determined from strength-duration 
measurements. 


(d) In the case of implant electrodes during chronic pacing, the variations in the parameters 
of the strength-duration curves for the same electrode-tissues-target element system could 
be so significant that its influence on the optimal pulse shape can’t be neglected (for 
example, when the distance between the electrode and the target element varies). 


(e) The present approach, related with the theory of matched filters in electric engineering, is 
simple and straightforward. But it is suitable only for the criterion of minimum dissipated 
electric energy per pulse, in ohmic tissues. In other cases, the more powerful methods of 
functional analysis must be used “”, 
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Introduction: The posed problem (to find an optimal cathodic pulse shape (PS) for chronic 
autonomous pacing at threshold a given excitable tissue) can be approached from two 
standpoints: minimizing electric energy dissipation in tissues or minimizing battery 
drainage. 


Materials and Methods: An dissipation functional (EF) is introduced and defined as the 
temporal convolution of an impulse response function (IRF) with a cathodic PS. Then a 
given cathodic PS is just threshold if and only if the value of the EF equals a threshold 
charge (Th. Ch.). Both, IRF and Th. Ch can be experimentally determined from measured 
strength-duration curves and some additional information. Electric energy dissipated in the 
tissues (ohmic resistance) and in the whole system (in series combination of ohmic 
resistance of electrode’s cable, tissues resistance and polarization resistance of electrode- 
tissue interface - given by an Abel’s operator) is described by the corresponding 
dissipation functional (DF) of the PS. 


Results: The variational problem of finding the cathodic PS that minimizes one or another 
DF meanwhile the EF remains at its threshold value (or fixed over it) is solved by methods 
of functional analysis. The obtained PS is a monotonic increasing function of time, in good 
accord with known experimental facts for cardiac pacing. Bioelectric properties of tissue’s 
target region -among other variables- determine optimal PS. 


Conclusions: Optimal cathodic PS is a property of the electrode-tissue system. For the 
same electrode-tissue system, it can vary during chronic pacing. But if wished, optimal PS 
can be found from measurements with rectangular pulses. 


Keywords: autonomous chronic pacing at threshold of an excitable tissue, optimal 
cathodic pulse shapes, autonomous chronic pacing of the heart, electric energy 
dissipation in tissues, battery drainage, ohmic resistance of electrode’s cable, tissues 
resistance, polarization resistance of electrode-tissue interface, dissipation functional, 
excitation functional, restricted variational problem, functional analysis methods 


(A) Introduction 


During the first stages in design and construction of implantable cardiac pacing systems, 
the problem of identifying pulse shapes and durations that minimize the charge drainage 
from the batteries was posed (Angelakos and Torres, 1964; Roy and Wehnert, 1971; 
Fiandra, 1985, chapter 7). 
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Now, this is not a pressing problem as it was before, but minimizing the energy dissipated 
in the tissues in conditions of acute and chronic stimulation continues to be a goal to 
achieve. 


The problem to be solved here, from a physical-mathematical point of view, is to find an 
optimum cathodic pulse shape for chronic autonomous pacing at threshold a given target 
element in an excitable tissue. 

The target element may be a nerve fiber or skeletal muscle fiber (Rattay, 1990) or an 
electrical syncytium such as myocardium or electrically coupled visceral smooth muscle. 
The response of the target element to the applied pulse is the all-or- none emergence of an 
action potential. 


In principle the optimization problem can be approached from two standpoints: minimizing 
electric energy dissipation in tissues or minimizing battery drainage. In this research, to 
simplify the mathematical analysis, we work with a linear model of the electrode-tissue 
system relating the voltage between an active electrode (cathode) and an indifferent one 
(the anode), with the electric current injected in the biological tissues through the 
stimulating system. 


From a phenomenological standpoint the excitation process of a given target element in an 
electrode-tissues system may be described by an excitation functional (Fiandra, 1985, 
chapter 8; Suarez-Antola, 1994). In this research, also to simplify the mathematical 
analysis, we use a linear approximation to the excitation functional. 


(B) Materials and Methods 


To solve the optimum problem, we introduce two functionals: an excitation functional y[i| 
and an energy dissipation functional g]i ib 


The excitation functional: For a cathodic make excitation with a pulse i(t) of durationt - 


the excitation functional gives a necessary and sufficient condition to identify a just- 
threshold cathodic pulse. 


This excitation condition, in the linear approximation, reduces to the following equality, 
where Q, represents a certain threshold charge, and G(t) is an impulse response function: 


vli]=(g.i) =] G¢,-1)-i@)-at=9, (1) 
0 
Here g(t)= Gl, -t) and (g,f) = {el: f(t): dt is an internal product in the space clot, | 
0 
of continuous real functions in the interval [0,1 |. 


The impulse response function G(t) is positive between t=0 and ¢=1, where f, is the so- 
called useful time up to the rheobase (see Figure 1) 
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Figure 1: The impulsive response function of a linear approximation to the excitation 
functional (Suarez-Antola, 1994). 


To apply Equation (1) as a starting point, we suppose that at the beginning of the current 
pulse the excitable membranes of the target element are at their rest state. 

Both, the impulse response function, and the threshold charge can be experimentally 
determined from measured strength-duration curves and some additional information. 


The energy dissipation functional: The energy dissipation functional gli can be obtained 
from the power dissipated by a cathodic pulse in the active electrode interface, in the 
electrode leads and, in the tissues, neglecting the counter-electrode effects. If V(¢)is the 
voltage between the input end of the active electrode lead and the indifferent electrode, then 


we put: gli]= i V(t)-i(t)-dt In the linear approximation (Suarez-Antola, 1994): 
0 


V (0) =R, -ilt)+ | K(e—u)-ilu) du (2) 


The ohmic resistance R, =R,+R, is the sum of an ohmic lead resistance R, and the 


tissues resistance R,. This last one, in the conditions of electric stimulation can be taken as 
ohmic also (Schaldach y Furman, 1975; Fiandra, 1985; Rattay, 1990). 


The polarization resistance of electrode-tissue interface is described, in the linear 


t 
approximation, by the convolution term | K (t - u): i(u)- du. 
0 


Then, for a cathodic make excitation with a pulse i (t) of duration ?, the energy dissipation 


functional gli | may be written, in the linear approximation: 
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ty 


gli|=R, Lesa fay| fst) aaa at 
0 0 oO 
In inner product notation we have: 


gli]= R, -(i,i)+ (Ki, i) (3) 
Here, by definition, the linear operator K verifies: 

Ki(t)= [KC —u)+i(u)- du (4) 
According to De Boer, R. and A. van Oosterom (1978) the kernel of this operator is given 


by K(t)=K, {é , where K,,7, and mare positive parameters (0<m<1). It is an 
0 

Abel’s operator. The properties of the interface between the active electrode and the tissues 

can vary during chronic stimulation (Schaldach, M. y S. Furman, 1975; Irmich, 1985; 

Rattay, 1990) so we can expect some variation in the above introduced parameters. 

The energy dissipation functional gli is positive definite because both dissipation terms, 

R, (i, i) and (Ki, i) , are strictly positive if the cathodic current pulse is not identically 


Zero. 


The optimum problem: To find a just-threshold cathodic pulse shape i, (t), continuous in 
the time interval lo, 1 ; | , such that: 


(1) y|i, |=, (just-threshold condition) while 


(II) gli, |< gli, +6 il for every other pulse shape i,+06i such that wii, +6 il= OQ, (so 
ij + 07 is any other just threshold pulse shape). 


Notice that, due to the linearity of the excitation functional, from w|i, |= QO. and 


y li, +61 | = Q, it follows that for every variation di : 


y|6i]=0 (5) 
From Equation (3) we derive: 


gli, + 5i]- gli. ]=2-R, -(ip,51) + (Ki, 51) + (i, K5i) + 0(53) (6) 


Here o(di)=R, -(5i,5i)+ (6 i,KO i) is always positive if Si is not identically zero. 
So, if we find a current just threshold pulse i, (t) such that for every perturbation di that 
verifies y[5i]=0, we have2-R,, -(i),51) + (Kip, 51) + (ig, K5i) =0 , then i,(t) will be an 


optimum pulse shape. 
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Now, the continuous operator K has an adjoint operator K* (Cotlar and Cignoli, 1974; 
Griffel, 1983). 


As consequence: (i),K5 i) = (Kiy,5i) (7) 
Taking (7) into account we derive: 2-R, - (i), 51) + (Ki,,5i) + (i), K5i) =2-R, -(Ai,,5i) 
Here, by definition, the operator A verifies: 

pes em | ae (8) 


Dp 


As usual, / is the identity operator. 


Considering that, from Equation (1) it follows wo i]= (g,6 i) and considering Equation 
(8), the optimum problem can be recast as follows: 
To find i, (t) such that: (iy, 6 i) =0 for every perturbation that verifies (g, 67) =0 


(C) Results and Conclusions 


If two linear functionals like (Ai, 53) and (g,di) have the same null set, the functions 


Ai, (t) and g(t) must be proportional (Cotlar and Cignoli, 1974): 


Ai,(t)=A- g(t) (9) 


Equation (9) is a convolutive Volterra’s integral equation of the second kind, so it always 
has a unique solution for every continuous function g(t) in lo, t oI (Moissiiewitsh, 1977). 


If 4-lis the inverse operator: i(t)=A-A g(t) (10) 


Now, as i, (¢) must be a just-threshold pulse: 


wlin]=(g,i,) = 9, (11) 
Substituting (10) in (11) we obtain: 
y= Ag (12) 
(g. A g) 


This is the optimum threshold current pulse that we were searching. 


It can be shown (De Boer and van Oosterom, 1978) that the effect of the polarization of the 
interface can be often neglected. In this case the shape of the optimum pulse is determined 
mainly by the function g(t)= Ge A t). 
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As will be studied elsewhere, known experimental results (Angelakos, E. and J. Torres, 
1964; Roy and Wehnert, 1971; Klafter, R. and L. Hrieben, 1976) agree with this prediction. 


The experiments show that the parameters of the strength-duration curves, and therefore the 
threshold charge and the impulse response function depend, besides the threshold properties 
of the membranes of the target element, of certain geometric and bioelectric parameters of 
the system formed by the electrodes, the tissues, and the target element within the 
biological tissues. For example, they depend on the distance between the cathode and the 
target element, and of the size and shape of the active electrode if it is located close enough 
to the target element. 

The ohmic term in the dissipated energy per pulse is proportional tok, =R.+ Rr. 


But R; usually depends on the size and shape of the cathode (we suppose monopolar 
stimulation) and it is very sensitive to conductivity variations in the tissues located near the 
active electrode, such as the ones related with inflammatory reactions to implant electrodes 
(Fiandra, 1985). 

Consequently, the optimal pulse shape for a given duration can be modified, in the same 
electrode-tissues-target element system, due to the variation of the parameters of the 
system. 


To end: 


(a) Optimal pulse shape is a property of the electrode-tissue-target element system. For the 
same electrode-tissue-target element system, it can vary during chronic stimulation of 
myocardium, nerve, or skeletal muscle. 


(b) In principle, optimal cathodic pulse shape could be found at each moment from 
measurements of the strength-duration curves and from a characterization of the 
polarization resistance of the active electrode interface, employing rectangular pulses of 
controlled electric current. For a given duration, and a just-threshold pulse, the optimal 
shape is given by equation (12). 


(c) In the case of implant electrodes during chronic pacing, the variations in the parameters 
of the strength-duration curves for the same electrode-tissues-target element system 
could be so significant that its influence on the optimal pulse shape can’t be neglected 
(for example, when the distance between the electrode and the target element varies). 


(d) A nonlinear analysis of the optimum cathodic pulse shape remains to be done. 
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Contributions to the study of optimal biphasic pulse shapes for 
functional electric stimulation: An analytical approach using the 


“excitation functional”. 


Roberto Suarez-Antola 


Abstract—An analytical approach to threshold problems 
in functional electric stimulation and pacing is proposed, 
framed in the concept of excitation functional. This 
functional can be applied to nerve, muscle, and 
myocardium stimulation by external electrodes. An 
optimal shape for a biphasic pulse is found, using the 
criteria of minimum energy dissipated in biological 
tissues and total charge compensation between the 
excitatory cathodic and the compensatory anodic phases. 
The method can be further developed and applied to 
other threshold problems in _ functional electric 
stimulation and pacing. 


I. INTRODUCTION 


ELECTRICAL stimulation of excitable biological tissues 
by working electrodes, located near the target elements 
in which action potentials must be produced, pose 
several interesting problems from an engineering 
standpoint. 

Reference [1] reviews the state of the art in relation 
with the design of efficacious and safe protocols for 
functional electric stimulation, both from an 
electrophysiological and clinical point of views, 
considering the electrochemical aspects of charge 
injection and the mechanism of damage to tissues and 
electrodes. Biphasic current pulses are needed. 
Usually, a cathodic half pulse of rectangular shape is 
followed by an anodic one (rectangular or 
exponentially decreasing), often after a short delay 
time. The charge compensation thus obtained at the 
electrode-electrolyte interface, allows pacing the target 
tissue without a significant damage to the biological 
tissues located near the working electrode or the 
counter electrode (necessary to close the circuit) [1], 


[2]. 


The purpose of the present paper is to contribute to the 
design of stimulation protocols that, besides being 
efficacious and safe, may be considered as optimal 
from an engineering point of view. Minimum charge 
drainage and transfer, minimum energy dissipation in 


biological tissues or in the whole working electrode- 
tissues-counter electrode system, as well as optimal 
durations of the cathodic and anodic components of the 
biphasic pulses and of the delay time between the 
cathodic and the anodic parts are some reasonable 
design aims that may be imposed for just threshold 
electric current pulses [3]. So, a threshold condition 
that allows us to classify biphasic pulses as below, 
above, or just at threshold is needed to discuss and 
solve optimization problems. After that, a safety factor 
must be introduced to apply the obtained results in 
clinical practice. 

Thresholds for different pulse shapes can be 
measured in true electrophysiological experiments in 
realistic settings. Complementary information may be 
produced by digital simulation using standard 
computing codes that allows us to consider the effect 
of an external working electrode in a target fibre or 
syncytium with varying geometries. But as our 
intention is to develop an analytical approach to 
optimization problems, we will use the excitation 
functional as main mathematical tool. In a previous 
work a simplified version of the excitation functional, 
(taking into account membrane accommodation 
processes but assuming that activation processes were 
instantaneously relaxed to equilibrium with membrane 
voltage), was used to study just threshold monophasic 
cathodic pulse shapes that minimize the energy 
dissipated in biological tissues [4]. As a more complex 
and realistic example of the proposed method, we 
discuss cathodic make excitation with biphasic pulses 
in which the anodic pulse achieves total charge 
compensation. We perform this task using a theoretical 
framework grounded in a more complete version of the 
functional, that also allows to study both one-phase 
and biphasic anode break excitation. The idea of an 
excitation functional is implicit in the early classical 
one and two factor models of excitation (due to 
Monnier, Hill, Rashevski, and generalized later by 
FitzHugh, although only for uniformly polarized 
membranes) [5] [6]. However, it seems that the first 
explicit proposal of the functional was done in a 
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contribution to the 1994-World Congress in Medical 
Physics and Biomedical Engineering [3]. 

Here we postulate the functional and its main 
properties. But it is possible to give to it a rigorous 
biophysical foundation grounded in equations of 
membrane dynamics and the modern tools of nonlinear 
science [7]-[9]. 


II. THE “EXCITATION FUNCTIONAL” 


A. Definition of the excitation functional 


In the following we consider in general the electric 
stimulation of nerve and muscular fibres, as well as 
myocardium and other electric syncytia. 

For each electrodes-tissues-target element system, 
we introduce a dimensionless impulsive response 
function G(t)that characterizes an equivalent charge 
transfer response to an applied electric current 
pulse i(t) . A formula for G(z) and its description will be 
given in the next section. 

From now on, cathodic pulses will be taken as 
positive functions of time, and anodic pulses as 
negative time functions. 

If the excitable membranes of the target element are 
in their rest state when ¢ <Oand the pulse begins 
whent = 0, let us construct the following equivalent 
charge: 


q.(t)= jolr—u)i(aa (1) 


Now, let us consider a pulse i(+)and let us take the 
maximum value of q,(t)for0<¢. This procedure 
makes a correspondence between each function i(t) 
and a real number: this defines a functional F|i]. 


There is a positive threshold value q,, such that if a 
pulse i(t)verifies Fli]< Yer it is below threshold (the 


action potential doesn’t emerge in the target tissue), if 
F\i|> q.rit is above threshold (an action potential 


emerges in the target tissue), and if F[i]= Jer it is a just 
threshold pulse. 

While the convolution that defines g,(t) is linear 
operation, the excitation functional Fli] is a nonlinear 


correspondence between functions and numbers, due to 
the appearance of the maximum of 4q,(t)in its 


definition. 
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B. The impulsive response functions. 
Then the impulsive response is: 


(~2# es") 


Pe a 
As - (A, + Ay) Ay — A (2) 


Here 4,,4,,A, are kinetic parameters of the electrode- 


tissues-target element system. (With the opposite sign 
convention for electric current pulses, —G(t) must be 


used instead of G(t)). The parameter A, is related with 


the activation processes in the excitable membranes, 
and it is at least one order of magnitude greater than 
both, 4,and4,. When the distance between the 


electrode and the target tissue is great enough in 
comparison with the space constant of the tissue, and 
the accommodation processes are fast enough, A, and 


A, are complex conjugate numbers. If not, they are real 
numbers and 4, always smaller thaniz,. Heaviside’s 
unit step function H(t) is equal to 0 if ¢<Qand to 1 
i>, 


Figure | shows a typical impulsive response function, 
when both parameters 4, and A, are real numbers. 


Function 


Impulsive Response 


time 


Fig. 1. Sketch of an Impulsive Response Function. 


G(t)is positive between Oandz¢,, and negative after. 
The growing of G(t), from G(0)=0 up to G(t,,)=1, is 
due to the activation processes in the excitable 


membranes of the target tissue. The negative part after 
t, 1S due to the accommodation processes in those 


membranes. 

Considering the difference in the time scales of 
activation processes, membrane time constant and 
accommodation processes, the following approximate 
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formulae are derived from (2): 


As 


= —lIn 
= As (A, +A) ©) 
1 A, 
ty, ¥—T lh = (4) 
(A, -A,) Ay 
tn QL, (5) 


If the activation process is instantaneous: 2, = 


Then formula (2) may be simplified: 


a= Ge A b*") = 4 le } a0 


When the accommodation processes are slow enough, 
A,is at least one order of magnitude smaller than 4, , 


and it is possible to make the approximation J, =0. 
Then: G(t)=e"”* A(t) (7) 


It is possible to express 7, as a well-defined function 


of the membrane parameters of the target element 
(including time and space constants), distance between 
the working electrode and the target tissue, and 
electrode’s radius [7]-[9]. 

Recent applications of (7) to the study of time 
constants for cathodic make stimulation of nerve and 
muscle fibers as well as electric syncytia like 
myocardium, may be found in [10] [11]. 


C. Application of the excitation functional to a 
cathodic pulse of constant amplitude: strength- 
duration curves for cathodic make excitation with 
monophasic current pulses 


Let us consider a rectangular pulse of durationt, and 
positive amplitude /, : 
i(t)=1, [A(e)- H(r-#,) (8) 


From (2) we obtain: 
t 


If <r, then: declt)= 1, Gu)du (9a) 

0 
d 

So: = = 1,.G(t) (9b) 

Is positive if0 <¢<z, and it is negative ift >¢,, (Fig.1). 
t 

If ¢. <¢ then: deelt)=1. [G(w)du (10a) 

t-t, 
d 
So: deel) _ (G4) -t.) (10b) 


Is positive ifG(t)>G(t-t.), it is zero in case of 
equality and negative if the inequality is reversed. 
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For fixed pulse amplitudes, the maximum value of 
dec(t) is attained whenr, =¢,. The minimum 


threshold amplitude (rheobase) I, », then verifies: 


t, 
TeRb a G(u)du — Ver 
0 


When ¢, <t,the maximum of q,., (t) is attained after 


(11) 


the end of the pulse, for an instant 7 such that (Fig2): 


Gl?)= Gr -1,) (12) 


Impulsive Response 
Function 


time 


Fig. 2 Sketch of the graphical meaning of equation 
(12). 


This is due to the increasing (between 0 andz,,) and 
decreasing (after t,, ) phases of G(+) 
If ty, =0 (instantaneous activation of the membrane) 


the increasing phase disappears and the maximum of 
Gog (t)is attained always at the end of the pulse. 


Solving (12) with (2), and considering the different 
time scales to make some approximations: 


The functione,(r,)decreases fromt,,to 0 when 
t. increases from: t, =0 to t, =t, 

The excitation functional gives the following analytical 
and threshold 


amplitudes /,, (strength-duration curve): 


Je, 
Tp (t.)= eur. , . (14) 


[Glu)du 


formula for pulses of duration?, 
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The limit of 7,.J,7(¢,)when 1,tends to Qis the so 
called cathodic limit threshold charge:Q.79 It 1s 


experimentally determined from the measured 
strength-duration curves. From (14) we find that 


dec 0..7.0 (15) 
From (11) we obtain this formula for the cathodic time 


constant of the 


(16) 


D. Thresholds for rectangular biphasic pulses with 
delay between the cathodic and the anodic parts 


system:t,5 = 


t, 
a [GW).du 
I 0 


c,Rb 


Let us consider a biphasic pulse i,(t) formed by the 


superposition of two rectangular component pulses. 
The first one,i,(t) cathodic, has positive amplitude 


I,and durationt,. The second one, i,(t)anodic, has 
negative amplitude /, and durationt,. There is a time 
delay 4, between them. The duration of the biphasic 
pulse is t, +tg +t, =t, 


The equivalent charge is the sum of two 
contributions: 
gelt)=deelt)+ dealt) 
(17) 
i. 
Ges (t)= eG —u)i.(u)du (18) 
0 
t 
Ge.q(t) = J G(t—u)i,(u)du (19) 
0 


ea(t)=Oif O<t<t,+tz, So that in this time interval 


gelt)=4eclt) 


After that, the negative anodic equivalent charge 
qe,(t)adds to the positive equivalent charge and tends 


to diminish q,(t) If ¢, <e.(¢.), the maximum of 


q(t) will be less than q, .(¢, + €,(t,))and to just reach 
qerthe amplitude of the cathodic pulse needs to be 


higher, raising the threshold. 


If t;>¢-,(¢.)this effect disappears: the anodic 


compensating pulse doesn’t interfere with the threshold 
of the excitatory cathodic part. 


Il]. BIPHASIC PULSES WITH CHARGE COMPENSATION 
AND LEAST ENERGY DISSIPATION IN BIOLOGICAL 
TISSUES 


Here we are going to apply the results obtained in 
Part II to an optimization problem: find the shape or 
shapes i,(t)=i.(r)+i,(t) of just threshold biphasic 
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pulses that produce a minimum of dissipated energy in 
ohmic biological tissues. In quantitative terms, we seek 
to make the quadratic functional: 


ly 
fi,’ (Cat 
0 


a minimum maintaining the excitation functional at 
its threshold value: 


(20) 


Flip ]=4er (21) 


ty t, tottg tte 
Note that [i,7(thdt=i?e(thde+ | 
0 


0 t.+ty 


i*a(t)dt (22) 


Besides the just threshold condition (21), we introduce 
two further restrictions. 
The first one is that ¢, should be greater than f,, so that 


we could expect that the anodic part doesn’t interfere 
with the threshold of the excitatory cathodic part, even 
if this last one doesn’t have a rectangular shape. In this 
case: F|i,|=F[i,] Then the threshold condition must 


be applied only to the cathodic part of the pulse. 


The second restriction is complete charge 
compensation: if i,.r(t)is the optimal just threshold 


cathodic shape for minimum energy dissipation, the 
anodic shape i,(t) should minimize energy dissipation 
while verifying: 


é 
Jing Qdt+ 
0 


to+ty +H, 


fi, (Qldt = 0 


t.+ly 


(23) 


Let us assume that g,.(¢), given by (18), takes the 


threshold value g,,; for t=t,+¢, In this case, 


applying § Cauchy-Buniakowski-Schwartz (CBS) 
inequality ([12]), we have: 
2 
t, 
qQ’er = [G¢. +é,—u)i,(u)du| < 
: (24) 


U+&, 


CJ G(w)duy(fie(u)du) 
& 0 


The minimum of dissipated energy is attained when the 
inequality becomes equality, and this happens for 
O0<¢<t, if and only if: 
i(t)=i.7(()=K.G(t, + €. -t) (25) 
The proportionality constant is calculated substituting 
(25) in the just threshold condition: ¢,.. (¢,+6,)= Ver 


We obtain: 
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Ve.T 
K.= 
t+& 


J G? (u)du 


Eo 


(26) 


to+tg tty 


Ji,(Qdt is 


to+ty 


Now, once i..r (t)is determined, Q, = 


fixed by (23). 
But if 1(t)=1for everyr, then applying a second time 


CBS inequality: 


2 


to+tg +t, t.+ty +H, 
| ae a ent ea { i, t .dt 
t.+ty t,+ty 


This can be equality, and the energy dissipated in the 
tissues during the anodic phase is a minimum, if and 
only if the anodic pulse is a rectangular one: 


i,(t)=K,{A(t-(t. +tg))-A(t-(t. +tg +t} (27) 


The restriction (23) gives: 


K,= © (28) 


a 
Thus, given the impulsive response function of the 
electrode-tissues-target element system, the durations 
of the cathodic and anodic phases, and a suitable lower 
bound for the interval of time between these phases, 
the shape of the just threshold biphasic pulse that 
minimizes the energy dissipated in the tissues was 
derived from the excitation functional and the well- 
known CBS inequality. 
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